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I. INTRODUCTION: DIFFUSION PROCESSES 

Among the elementary processes that underly natu- 
ral phenomena, diffusion is certainly one of the most 
ubiquitous. In an ensemble of moving elements -atoms, 
molecules, chemicals, cells, or animals- each element usu- 
ally performs, at a mesoscopic description level, a random 
path with sudden changes of direction and velocity. As 
a result of this highly irregular individual motion, which 
is microscopically driven by the interaction of the ele- 
ments with the medium, and of the elements with each 
other, the ensemble spreads out. At a macroscopic level, 
this collective behavior is -in contrast with the individual 
microscopic motion- extremely regular, and follows very 
well defined, deterministic dynamical laws. It is precisely 
this smooth macroscopic spreading of an ensemble of ran- 
domly moving elements that we associate with diffusion. 

One of the first systematic observations of diffusion 
was made by the botanist Robert Brown in 1828. He 
noticed that pollen particles dispersed in water exhibit 
a very irregular, swarming motion. In 1905, Einstein 
conjectured that this "Brownian motion" is due to the 
interaction of pollen with the water molecules and, in 
fact, proved that microscopic particles suspended in a 
liquid "perform movements of such magnitude that they 
can be easily observed in a microscope, on account of the 
molecular motions of heat" [|l[. Since then, Brownian 
motion is used as a synonym of diffusion. 

A very suitable and very useful mathematical model 
for Brownian motion is provided by random walks [^. 
In its simplest version, a random walker is a point par- 
ticle that moves on a line at discrete time steps At. At 
each step, the walker chooses to jump to the left or to 
the right with equal probability, and then moves a fixed 
distance x. This stochastic process can be readily gener- 
alized, firstly, by allowing the walker to move in a many- 
dimensional space. In addition, time can be made contin- 
uous by associating a random duration with each jump 
or by introducing random waiting times between jumps. 
Finally, the length of each jump can be also chosen at 
random from a continuous set with a prescribed proba- 
bility distribution. 



Being a stochastic process, a random walk admits a 
probabilistic description in terms of probability distribu- 
tions for the relevant quantities |3|. In particular, one 
is interested at studying the probability of finding the 
walker in a certain neighborhood dv of point r -in gen- 
eral, in a d-dimensional space- at time t, P{r, t) dr. Note 
that, besides its interpretation as a probability distribu- 
tion, P{r, t) can be related to the density in an ensemble 
of noninteracting identical random walkers. In fact, if the 
ensemble contains N walkers, n(r, t) = NP{r, t) stands 
for the space density of walkers. 

Suppose that the random walk is defined in continuous 
time and space, with a waiting time probability distribu- 
tion iP{t) and such that the probability that the walker 
jumps from any point r to r -|- x is p(x) dx. The normal- 
ization of probabilities imposes 



'0(t) dr — 1, 



J p{x) dx ~ 1. 



(1) 



If, moreover, ipir) and p(x) satisfy 

(r) = y T ipir) dr < 00, {x^) — j x^p(x) dx 



< 00, 



(2) 



[x = |x|) it can be proven that P(r, t) obeys the diffusion 
equation Q 



dP 

at 



= DVf P, 



(3) 



where D oc (a;^)/(r) is the diffusion constant, or diffusiv- 
ity. This equation has to be solved for a given initial con- 
dition P(r, 0) with suitable boundary constraints. The 
density n(r, t) of an ensemble of noninteracting diffusing 
particles obeys the same equation. 

A typical solution to the diffusion equation describes 
a density profile that, as time elapses, is smoothed out 
and broadens. In fact, it can be straightforwardly shown 
from the general solution that the width of the spatial 
distribution grows with time as 



j r^P{r,t) dr ^ 2dDt, 



(4) 
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where d is the space dimension. Correspondingly, it can 
be shown that the mean square distance between the 
present position and the initial position of a random walk 
that satisfies Eq. (H) is proportional to time. This pro- 
portionality between mean square displacement and time 
is the fingerprint of diffusion, as it can be used experimen- 
tally, numerically, and theoretically to detect this kind of 
transport mechanism in a given natural process. 

Though, being a form of transport, diffusion is inher- 
ently a nonequilibrium process, the large-time asymp- 
totic dynamics of an ensemble of diffusing particles can 
be described in the frame of equilibrium statistical me- 
chanics. In fact, it is expected that, for very large times, 
the system reaches a state of thermodynamical equilib- 
rium with the medium -and between the particles, if they 
interact. In such state, the diffusing particles and the 
medium participate of a balanced interchange of momen- 
tum and energy, which mantains the particles in their 
characteristic irregular motion. Once this situation is 
reached, a connection between the parameters that char- 
acterize thermodynamical equilibrium and particle dy- 
namics should exist. Einstein investigated this problem 
in 1905 ||l|, and concluded that diffusivity and tempera- 
ture are proportional: 

D = fikeT. (5) 

Here ks is Boltzmann constant, and /i is the mobility. 
The mobility is defined as the inverse of the friction co- 
efficent, in the present case, of the diffusing particles in 
the medium . The Einstein relation, Eq. (|^) , provides 
thus the expected connection between diffusion and ther- 
modynamical equilibrium. 

Despite the omnipresence of diffusion as a transport 
mechanism in natural processes, a different kind of trans- 
port underlies a selected -but ever growing- class of sys- 
tems. Due to various motivations most of these sys- 
tems have recently attracted very much attention. They 
range from turbulent fluids, to chaotic dynamical sys- 
tems, to genetic codes (see next section). In these sys- 
tems, anomalous diffusion -a mechanism closely related 
to normal diffusion, but with some qualitatively different 
properties- drives transport processes [||. Over the last 
few years, it became more and more clear that anomalous 
diffusion can be made naturally compatible with equilib- 
rium thermodynamics if the Boltzmann-Gibbs formula- 
tion of thermodynamics is replaced by Tsallis'. This com- 
patibility generalizes then the connection between normal 
diffusion and the usual formulation of thermodynamics. 
The main aim of this paper is to review this general- 
ization, commenting on some additional related topics 
brought to light in recent work. 

II. ANOMALOUS DIFFUSION AND LEVY 
FLIGHTS 

Any transport mechanism which, like diffusion, be- 
haves at mesoscopic level as an isotropic random pro- 



cess, but which violates Eq. is generally refered to 
as anomalous diffusion. More specifically, most of the lit- 
erature on anomalous diffusion has been devoted to pro- 
cesses where the mean square displacement (r^) varies 
with time as 

(r2) oc t^/', (6) 

where z 2) is the dynamic exponent, or random walk 
fractal dimension, of the transport process. Normal dif- 
fusion corresponds to z = 2. For z > 2, the growth 
rate of the mean square displacement is smaller than in 
normal diffusion, and transport is consequently said to 
be subdiffusive. On the other hand, for z < 2 the mean 
square displacement grows relatively faster and transport 
is thus superdiffusive. In the following, the attention will 
be mainly focused on this latter case. 

A. Anomalous diffusion in Nature 

As advanced above, anomalous diffusion occurs in a 
wide class of natural systems and processes. In the realm 
of physics, a paradigmatic example is given by particle 
transport in disordered media. Consider the motion of 
particles in a medium containing impurities, defects, or 
some kind of intrinsic disorder, such as in amorphous ma- 
terials. Examples are disordered lattices, porous media, 
and dopped conductors and semiconductors. In these 
heterogeneous substrates, particles are driven by highly 
irregular forces, which determine a complex variation of 
the local transport coefficients. This heterogeneity can in 
fact induce anomalous diffusion. For instance, it has been 
experimentally shown that in quasi-one-dimensional ionic 
conductors such as hoUandite (Ki.54Mgo.77Ti7.230i6), 
where transport is very sensitive to the presence of im- 
purities, the dynamic exponent is given by 

+ i (7) 

with 9 proportional to the temperature ||^. 

A reasonable model for transport in heterogeneous 
media is provided by a random walk in a lattice with 
quenched disorder § . This disorder applies to the depth 
of the potential wells at each lattice site, and to the 
potential barriers between sites. Randomness in these 
parameters induces a distribution for the time that the 
walker spends at each site before hoping to a neighbor. 
Generally, this waiting time distribution, iP{t), behaves 
as iI^{t) ~ r~^~^ for r — > oo. For /i > 1, the first relation 
in Eq. holds, and normal diffusion is observed. On 
the other hand, for fj, < 1 the mean waiting time diverges 
and diffusion is anomalous. The corresponding dynamic 
exponent for < /i < 1 is [Q 

z = 2/fi (d > 2), z ^2-d + d/fi {d<2). (8) 

A most important instance of anomalous diffusion in 
physics occurs in turbulent flows. In fully developed tur- 
bulence, fluid particles exhibit very irregular motion over 



2 



a wide range of space and time scales. Based on empiri- 
cal motivations, L. Richardson proposed in 1926 that the 
probability P{R, t) that two fluid particles initially close 
to one another have a separation R at time t obeys the 
equation ||] 



dP 
'dt 



d_ 

OR 



, dP 



(9) 



Comparing with (||), it is clear that the Richardson equa- 
tion is a diffusion equation with space-dependent diffu- 
sivity. Its solution immediately implies {R?) oc i'^, indi- 
cating that the relative motion of particles in fully de- 
veloped turbulence corresponds to anomalous diffusion 
with a dynamic exponent z ~ 2/3, which is well into the 
superdiffusive regime. 

Richardson's law has to be modified to take into ac- 
count the fact that the vorticity field in a turbulent flow 
is intermittent This means that vorticity -and, in 
particular, turbulent activity and dissipation- is concen- 
trated on a relatively small volume in the whole sys- 
tem, which happens to be a fractal set. Experiments 
pof suggest that the fractal dimension of this set is 
df = 2.8 ± 0.05. Incorporating this correction, it can 
be shown [|ll that (i?^) oc ^12/(1+^/)^ gj. 

1 + ^/ 
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0.63. 



(10) 



A somewhat more abstract form of anomalous diffu- 
sion is present in the evolution of chaotic Hamiltonian 
dynamical systems in phase space. Hamiltonian systems 
are characterized by volume conservation in phase space, 
as stated by the Liouville theorem. The domain occupied 
by a given set of initial conditions in phase space can be 
strongly distorted under the effect of evolution, but its 
volume remains constant. Mechanical processes that pre- 
serve energy are instances of Hamiltonian systems, but 
a huge host of systems -both continuous and discrete 
in time- is known to belong to the same class. Since a 
single Hamiltonian system can exhibit both regular and 
chaotic evolution by simply changing the initial condi- 
tion, the dynamical geometry of its phase space is usually 
extremely intrincate. Zones of nested regular trajectories 
which alternate with chaotic regions are typically found 
at many scales, displaying selfsimilar structures. In the 
bulk of chaotic regions trajectories are extremely irreg- 
ular and resemble random paths. On the other hand, 
when approaching the boundary with a regular region, 
the same trajectory can temporarily become much sim- 
pler and smoother. Consequently, as it evolves in phase 
space along a chaotic orbit, a Hamiltonian system al- 
ternates intermittently between zones of highly complex 
behavior and a regime of almost regular dynamics. Glob- 
ally, this motion can be thought of as a stochastic process, 
and turns out to have the same statistical properties as 
anomalous diffusion. 

A case studied in detail in the literature is the so-called 
Q-flow . It is defined as a three dimensional velocity 
field 



X — d'i' /dy + esinz 
y= d'^/dx + ecosz 
z = ^ 



(11) 



with 



*(.T, y) = ^ cos [x cos(27rj/A:) + y sin(27rj/A;)] . (12) 



Here e is a parameter and k is an integer that determines 
the symmetry of the flow. The solution to Eqs. (|Tl| ) 
is a complex trajectory that wanders in an infinite con- 
nected net of channels of width of order e, inside which 
the trajectory looks like a random contour ||l^. Numer- 
ical measurements of the statistical properties of these 
diffusion-like trajectories show that the dynamic expo- 
nent z to be associated with them fluctuates strongly as 
a function of e ||l^. For fc = 6 and 0.8 < e < 1.8, z varies 
in the interval 



1 < z < 2, 



(13) 



making apparent that the motion is supperdiffusive, as 
in turbulence. Anomalous diffusion has also been ob- 
served in dissipative (non-Hamiltonian) dynamical sys- 
tems, both in simulations and in experiments. For exam- 
ple, a dynamic exponent z « 1.2 has been measured in 
the Taylor-Couette flow Q . 

As stated in the Introduction, anomalous diffusion is 
not restricted to physical systems. This kind of trans- 
port has in fact been detected to underly several biolog- 
ical processes. Some sectors of genomic DNA sequences, 
for instance, are known to exhibit statistical properties 
analogous to anomalous diffusion. To stress this corre- 
spondence, "DNA walks" have been deflned DNA, 
which codes genetic information, is a large molecule in 
the form of a chain of nucleotides. Each nucleotide con- 
tains either a purine or a pyramidine base. Two purines 
-the adenine (A) and the guanine (G)- and two pyra- 
midines -the cytosine (C) and the thymine (T)- are in 
turn present in the DNA chain. Therefore, the informa- 
tion code in DNA is a symbolic chain of four letters: A, 
C, G and T. Amazingly, within DNA only a small por- 
tion does code information for protein building (3% in 
the human genome) whereas other zones are noncoding, 
their specific role being unknown. A one-dimensional 
DNA walk is constructed by sequentially running over 
the chain of nucleotids. Each time a purine is found the 
walker jumps rightwards, whereas when a pyramidine is 
found the jump occurs leftwards. For instance, 

• • • ACGCTGAGTG • • • ^ ••• + - + -- + + + - + ■•• 

where -I- stands for jumps to the right and — stands for 
jumps to the left. In this DNA walk, systematic devi- 
ations from normal diffusive behavior derive from long- 
range correlations in the nucleotide sequence. It has been 
found that the DNA walk is in fact statistically identical 
to normal diffusion in the zones of the genomic chain that 
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code information. On the other hand, in noncoding se- 
quences the DNA walk is analogous to superdiffusion. In 
the human beta-globin chromosomal region the dynamic 
exponent is z « 1.4 ||l6|] . 

A less involved instance of anomalous diffusion in bi- 
ology appears in the flight patterns of certain birds. In 
particular, it has been found JTzt that in the foraging 
behavior of the wandering albatross {Diomedea exulans) 
the flight-time intervals exhibit a power-law distribu- 
tion. This results in a anomalous diffusive-like motion 
which, according to field measurements on the Bird Is- 
land, South Georgia, has a dynamic exponent z k, 1.2. 
The kind of flight patterns observed in these seabirds is 
supposed to reflect a complex structure in the underly- 
ing ecosystem, especially, in the spatial distribution of 
the exploited environment. It has been suggested that 
albatrosses specialize in long journeys of random forag- 
ing, searching for patchily and unpredictably dispersed 
prey over several million square kilometers. As discussed 
in the next section, anomalous diffusion is inherently re- 
lated with fractal geometry, scale-invariance, and self- 
similarity, which in the present example seem to drive 
the predator-prey dynamics. 

Of course, the previous collection of examples of 
anomalous diffusion in Nature is not at all exhaustive, 
but only pretends to give a hint on the variety of systems 
driven by this kind of transport. For a more detailed ac- 
count, the reader is refered to the review by Bouchaud 
and Georges ||^ . This review is also an excellent reference 
to the mathematical treatment of anomalous diffusion, 
which is briefly introduced in the following. 



B. Random- walk models of anomalous diffusion 

In view of the efficacy of random walks in modeling 
normal diffusion at a mesoscopic level, it is desirable to 
find a similar stochastic model describing anom alous dif- 
fusion. It has already been mentioned in Section II A that 
introducing a waiting time distribution which, for large 



waiting times, behaves as ipiT) t~ with < 1, pro- 
duces the anomalous dynamic exponent given in Eq. (|^). 
In this case z > 2 and, therefore, transport is subdiffu- 
sive. Note that the source of anomaly in these random 
walks is the divergency of the average waiting time (r) . 
For /i < 1, the long-tailed distribution ipiT) allows for 
very long waiting times with relatively high probability, 
violating thus the first relation in Eq. (||). These long 
waiting times produce an overall reduction of efficiency 
in the transport mechanism with respect to normal dif- 
fusion, and leads to subdiffusivc behavior. 

Taking into account the previous argument, in can be 
expected that the violation of the second of relations (|^) 
will lead, on the other hand, to superdiffusion. In fact, 
having a divergent mean square displacement {x^) re- 
quires the jump probability distribution p(x) to have a 
long tail for large a;, in a sense to be made precise im- 



mediately. This long-tailed distribution would produce, 
with relatively high probability, very long jumps. Glob- 
ally, this transport mechanism should result more effi- 
cient than normal diffusion, and superdiffusive behavior 
is expected. 

It can be easily shown that, in c?-dimensional space, 
the mean square displacement diverges if 



p(x) 



1 



(14) 



for large x, with 7 < 2. Note that, for such distribu- 
tion, normalization requires 7 > 0. It is therefore to be 
expected that a random walk with a jump distribution 
p(x) ^ x~'^~'^ with < 7 < 2 does not model normal 
diffusion, but some kind of superdiffusive motion. Of 
course, a large score of functions satisfy Eq. (p^), and 
are thus candidates to play the role of a jump distri- 
bution for a superdiffusive random walk. Among them. 
Levy distributions have been studied in detail. 

Levy distributions are defined through their 

Fourier transform, which reads 



p(k) 



j exp(jk • x) p(x) dx = cxp(— &fc^). 



(15) 



where & is a positive constant, and k = |k|. Although the 
antitransform p(x) has no analytical expression, it can 
be shown that it satisfies Eq. ([14]). Moreover, if 7 < 2 
the positivity of p(x) is insured. The relatively simple 
form of this jump distribution in the Fourier representa- 
tion makes it an ideal tool for analytical manipulation. 
However, the main interest of Levy functions in the math- 
ematical theory of distributions comes from the fact that 
they are stable. Essentially, this means that two Levy 
functions with the same Levy exponent 7 produce, upon 
convolution, a third Levy function with the same expo- 
nent. This can be readily proven in the Fourier represen- 
tation, where the convolution transforms into ordinary 
product: 



pi(k)p2(k) = exp(— 6ifc''') exp(-&2fc''') 

= exp[-(6i + 1)2)^^] = P3(k). 



(16) 



Levy functions are not the only stable distributions, 
the Gaussian p(x) cx exp(— x^) being probably the best- 
known example. The (one-dimensional) Cauchy distri- 
bution, p{x) oc (1 -I- a;^)~^, is another instance. Stable 
distributions play a fundamental role in probability the- 
ory since according to the central limit theorem -which 
is usually stated for the Gaussian function- the addi- 
tion of random variables tends to one such distribution. 
In particular, as P. Levy demonstrated through hisgcn- 
eralization of the Gaussian central limit theorem ||l8|,|l9| , 
adding random variables with a power-law distribution as 
in ( |l4| ) -whose second moment (a:^) diverges for 7 < 2- 
leads asymptotically to a Levy distribution. 

Another important property of Levy distributions, 
which is reflected in its power-law large- a; asymptotic be- 
havior, Eq. (14), is the absence of characteristic length 
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scales. This implies that random walks with Levy jump 
distributions have self-similar properties. In particular, 
it can be shown that the set of points visited by this kind 
of random walk is a fractal of dimension 7 ||2^ . As a con- 
sequence, these distributions ubiquitous in the realm of 
self-similarity geometry -the geometry of fractals. 

A discrete-time random walk whose jump distribution 
is given by a Levy function as in Eq. ( |l5| ) is called a Levy 
flight. It has been suggested |^ that the Fourier trans- 
form P(k, t) of the probability distribution of finding the 
walker at a given point at time t satisfies the evolution 
equation 



dP 
'dt 



(k,t) = -D-,k'^Pik,t). 



(17) 



This equation generalizes, in the Fourier representation, 
the diffusion equation (H). A straightforward dimen- 
sionality analysis shows that the anomalous diffusivity 
is D-y oc b/At, where At is the time step of the random 
walk. In free space, equation ( pT| ) can be readily solved: 

P(k, t) = P(k, 0) eyipi- D^k'^t). (18) 

For a delta-like initial condition, P(r,0) = S{r), one has 
P(k, 0) = 1 and, thus, 



P{r,t) = (27r 



exp(-ik • r - D^k'^t) dk. (19) 



This function remains unchanged, except for a constant 
factor, if both space and time are conveniently rescaled, 
P{a^^''r, at) = a~'^^''P{r, t). Using this scale invariance, 
one can write 



P(r,t) = t-''^m{r/t''^), 

where 11 is a function of a single variable 
this implies 

(r2) oc ^2/7^ 



(20) 
In turn, 

(21) 



for < 7 < 2. This result, which has been here derived 
for a delta-like initial distribution, can be generalized by 
simple superposition to more general initial conditions. 
It shows that a Levy flight with < 7 < 2 represents 
superdiffusion with a dynamic exponent 2 = 7. On the 
other hand, for power-law jump distributions with 7 > 2 
the dynamic exponent corresponds to normal diffusion, 
z = 2 (Figure |l|). 

The fact that in Levy flights the mean square displace- 
ment {x^) of a single step diverges, implies -in contrast 
with Eq. (pT|)- that the mean square displacement of 
the walker after a certain time is, on the average over 
infinitely many realizations, also infinite. The arguments 
used to derive Eq. (|2|) are therefore of limited valid- 
ity (2~~ 
result 



and have to be taken cum grano salts. The 
1) is expected to be valid for finite times, i.e. 
in a certain portion of the whole random walk, and on 
averages over a finite number of trajectories. The same 
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FIG. 1. The diffusion dynamic exponent z as a function of 
the Levy exponent 7. 

result would be valid during a certain time if the jump 
distribution p(x) is a Levy function in some (large) range 
of values of x, but has a cutoff for sufficently large x [g3| . 
In spite of this drawback, Levy fiights provide a very 
powerful tool for modeling superdiffusion because of the 
mathematical properties of the Levy distributions, sum- 
marized above. They are thus a very satisfactory starting 
point as a model for studying the statistical mechanics of 
superdiffusive transport, generalizing the results outlined 
in the Introduction for normal diffusion. 

Before passing to the discussion of superdiffusion in 
a statistical-mechanical frame, a comment is in order 
on the numerical simulation of superdiffusive random 
walks. Due to the cumbersome properties of Levy distri- 
butions in real space pO| ] , it is not convenient -in numer- 
ical calculations- to work directly with these functions. 
Rather, power-law distributions with the same asymp- 
totic properties as Levy's, Eq. (|lj), are used. For in- 
stance, one can take 

with M a normalization constant. The Fourier transform 
of these distributions behaves precisely like a Levy dis- 
tribution for small k, p(k) « 1 — bk'' . Also, they show 
the same scale invariance for large x, wich leads to self- 
similar properties in the associated random walks. In 
Figure the first 10^ points visited by a random walk in 
two dimensions, with the jump distribution given in (j2^) 
and 7 — 1.5, are shown. Note the clustered, fractal-like 
structure of this set of points. 



III. MAXIMUM-ENTROPY FORMALISM FOR 
ANOMALOUS DIFFUSION 

Entropy plays a central role in the foundations of equi- 
librium and nonequilibrium statistical mechanics. It is 
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FIG. 2. The first IC* points visited by a two-dimensional 
random walk generated by a power-law jump distribution 
with exponent 7 = 1.5, starting at the center of the main 
frame. The amplification illustrates the self-similar proper- 
ties of this process. 



well known from the work by L. Boltzmann and others 
that entropy provides a natural link between nonequilib- 
rium processes and their asymptotic states of thermody- 
namical equilibrium. In addition, the whole theory of 
equilibrium statistical mechanics can be derived from a 
variational formalism for the entropy, as follows. Define 
the entropy S* as a functional of the probability distribu- 
tion Pi over the states z of a given system, 



S[P] 



(23) 



where fcs is Boltzmann constant. Find then the values of 
Pi that maximize S[p], taking into account the normaliza- 
tion constraint, J2iPi — 1> ^^'^ required by particu- 
lar conditions of the system under study- any additional 
constraint on pi. The value oipi resulting from this max- 
imization procedure gives the probability of finding the 
system in state i when thermodynamical equilibrium has 
been reached. For instance, introducing the canonical 
constraint ^ eiPi = E, where is the energy of state i 
and E is the thermodynamical energy, the maximization 
of entropy produces the well-known Boltzmann distribu- 
tion Pi (X exp(-/3ei) ||]. 



A. Traditional formalism 

As a starting point for including normal diffusion in 
the frame of equilibrium statistical mechanics, the proce- 
dure of entropy maximization has been applied to obtain 
the jump probability distribution p(x) in a discrete-time 
random walk |^^. In this case, entropy is defined as a 
straightforward generalization of (p3). 



S[p] = -kB yp(x)ln[a'^p(x)] dx. 



(24) 



Here cr is a characteristic length, whose meaning will be- 
come clear immediately. The distribution p(x), in fact, 
has units of length to the power —d. The maximization of 
S[p] is carried out taking into account the normalization 
of p(x). 



p(x) c?x = 1, 



(25) 



and imposing the additional constraint 

J x^p{yi) dx = a^d, (26) 

which is inspired in the second relation of Eq. (^) . Except 
for a dimensionality factor, is thus the mean square 
displacement associated with p(x). 

Under these conditions, the maximization of entropy 
yields 



p(x) = {2t:(7^)''^/^ exp{-x^ /2a 



(27) 



namely, a Gaussian jump distribution. Since, in view of 
the constraint (p^), the mean square displacement associ- 
ated with p(x) is finite, the maximum-entropy formalism 
applied as above to the jump distribution of a random 
walk describes normal diffusion. 

The question on whether anomalous diffusion can be 
derived from a variational formalism for the entropy 
arises now quite naturally. MontroU and Shlesinger p^ ] 
have shown that this is in fact possible, but requires re- 
placing the constraint ( p6| ) by a more complex condition 
on p(x). In particular. Levy flights, Eq. ( [l^), a re ob- 
tained from the maximization of the entropy (pj) if the 
jump distribution satisfies, along with normalization. 



J In {2n)-'^ J exp{-ik-x^bk'^) dk 



p(x) dx 



= constant. 



(28) 



This is however a quite unsatisfactory answer to the 
above question. Indeed, besides its complexity, the con- 
straint ( p8| ) is anything but a natural condition to impose 
to the jump distribution. In MontroU and Shlesinger's 
words, "it is difficult to imagine that anyone in an a pri- 
ori manner would introduce" such a condition for max- 
imizing the entropy with respect to p(x). This remark 
would at once exclude Levy flights -and anomalous diffu- 
sion with them- from the frame of the maximum-entropy 
formalism and, therefore, from a natural connection with 
equilibrium statistics. 

In Ref. |2|] , a different approach has been proposed to 
tackle the problem of deriving anomalous diffusion from 
the maximization of entropy. Since replacing the con- 
straint on the jump distribution implies imposing uncon- 
ventional, forced conditions on p(x), a possible way out 
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is to replace the form of the entropy instead. In par- 
ticular, it has been found that the form of the entropy 
proposed by Tsallis produces, upon maximization 

with the constraints prescribed by this generahzed the- 
ory, power-law jump distributions with the asymptotic 
behavior given in (|l4|). As described in the following, 
random-walk models of anomalous diffusion find thus a 
natural statistical-mechanical basis in Tsallis' theory. 



B. Generalized formalism 

Inspired in the theory of multifractals, Tsallis |^ pro- 
posed to generalize the traditional Boltzmann-Gibbs sta- 
tistical mechanics by introducing new forms for the en- 
tropy and for the constraints to be applied in the maxi- 
mization procedure. For a system whose i-th state is oc- 
cupied with probability pi , the generalized entropy reads 



l-Q 



(29) 



where g is a real parameter. For the canonical ensemble, 
where the energy of the i-th state is and the average 
energy is Eg , the generalized constraint to be imposed to 
Pi, along with probability normalization, is 



E„ 



(30) 



In this generalized formalism, in fact, the average of any 
observable O is defined as {0)q = J2i OiPl- This average 
is usually refered to as the g-expectation value of O [|8| . 

The generalized statistical-mechanical formalism based 
on Eqs. (j^) and (|30| ) has some remarkable properties. 
First of all, it reduces to the traditional Boltzmann-Gibbs 
formulation in the limit g ^ 1. In fact, Eq. ( p3| ) is re- 
covered from (^) in that limit except for the factor fcs, 
which has here been conventionally put equal to unity. 
The canonical constraint (|3^) reduces in turn to the tra- 
ditional definition of mean energy. The new formalism 
preserves the full Legendre-transformation structure of 
thermodynamics for all q p7[ |, leaving invariant in form 
the main results of statistical thermodynamics, such as 
the Ehrenfest theorem, the H-theorem, the von Neumann 
equation, the Bogolyubov inequality, and the Onsager 
reciprocity theorem ||2^]. Its seems to be particularly 
useful in dealing with systems involving long-range cor- 
relations and non-extensivity, as the formalism itself is 
non-extensive for q ^ 1. The Tsallis exponent q has 
thus been interpreted as a measure of non-extensivity. 
Since its introduction a decade ago |^ Tsallis statistics 
has found successful applications to a large class of prob- 
lems of high interest, ranging from gravitational systems, 
to turbulent flows, to optimization algorithms. Many of 
these applications are described in detail in other papers 
of the present issue, and are therefore no longer discussed 
here. 



In order to apply Tsallis statistics to discrete-time ran- 
dom walks in the spirit outlined in the previous section, 
Eqs. ( p^ ) and (|2^) have to be generalized according to 
( ^ ) and (|30|), respectively. As a function of the jump 
pro bability, the generalized entropy can be written as 

Mull 



Sq[p] 



^jl-a-'l [aVx)]'dx|, (31) 



whereas the canonical constraint transforms into a con- 
dition on the g-expectation value of x^: 



[ap{x)]Ux^a^ 



(32) 



Here, a preserves its identification as a typical length 
associated with the jump probability. However, for q ^ 1, 
does not coincide with the mean square length of the 
jumps. For simplicity, the dimensionality factor in the 
right-hand side of Eq. ( p6| ) has now been absorved by a. 

It is shown in the following that the maximization of 
S[p\ as defined in (|3l|) with the constraints (^ and ( |3^ ) 
-which, as in the case of the traditional formalism, is car- 
ried out by the standard method of Lagrange multipliers 
|p6| , p7| - produces a power-law jump distribution. The ex- 
ponent of the power-law depends on the space dimension 
and on the Tsallis exponent q. This jump distribution is 
not a Levy distribution like (|l5|) , but has the same type 
of asymptotic behavior, Eq. (|14|). For suitable values of 
d and q, this form of p(x) will therefore define a random 
walk with anomalous properties. 

For the sake of clarity, the results in one dimension are 
shown first [ p3|j29| . The jump distribution resulting from 
the maximization procedure is, in this case, 

p{x)^Z-'[l-f3{l-q)x^f''-'^ (33) 
where the partition function Zq is given by 



Z„ 



-|-oc 



[1 - /3(1 - q)x' 



11/(1-9) 



dx. 



(34) 



The positive constant f3 is one of the Lagrange multi- 
pliers, which can be expressed as a function of a using 
the constraint (^2|), as shown below. In the generalized 
formulation of statistical mechanics /3 is related to the 
temperature T in the standard form, /3 oc 1/T. 

It turns out from Eq. (33) that the normalization of 
p{x) can be satisfied for g < 3 only. The Tsallis expo- 
nent for random walks is therefore restricted to the inter- 
val (— oo, 3). Explicitly calculating the partition function 
yields 



p{x) 



Piq - 1) 



(^) 



[l + (3{q~l)x^ 



-1/(9-1) 



(35) 
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FIG. 3. Jump probability distribution derived from Tsallis 
statistics, for some values of the Tsallis exponent q. For q = \ 
the standard Gaussian is obtained. Note the cut-off for g < 1, 
and the long-tailed power- law distributions for g > 1. 



for 1 < 5 < 3, and 



p{x) 



/3(1 - q)x^ 



ifx2 <l/;3(l-g), 



otherwise 



(36) 

for q < 1. For g ^ 1, of course, the Gaussian ( |27| ) is 
reobtained. Note that the Tsallis exponent q is related 
to the exponent 7 in Eq. ( jl^ ) according to 



7 = 



3-q 



3 + 7 
1+7' 



(37) 



Figure || shows the profile of p{x) for several values of q. 

Regarding anomalous diffusion, thus, it is clear that 
the case q < 1 is irrelevant. In fact, for such val- 
ues of the Tsallis exponent p{x) exhibits a cut-off at 
I a; I = 1/ — q), and vanishes for larger |a;|, as shown 
by Eq. (pw. This implies at once that the mean square 
displacement associated with p{x) is finite and the result- 
ing random walk corresponds to normal diffusion. The 
attention is consequently focused in the following on the 
case 1 < q < 3, Eq. (^5|). In this case, the mean square 
displacement is 



[/3(5-3g)]-i ifg<5/3. 



00 



iiq> 5/3. 



(38) 



Therefore, for 1 < q < 5/3 ~ 1.67 the mean square dis- 
placement is still finite, and the random walk corresponds 
to normal diffusion. On the other hand, anomalous su- 
perdiffusion is obtained for 5/3 < 9 < 3. 



It is interesting to calculate now the g-expectation 
value of x^ which, in the frame of Tsallis statistics, re- 
places -as an average quantity- the mean square dis- 
placement of the standard formulation. According to the 
constraint (j3^) imposed to p{x) in the maximization of 
entropy, this ^-expectation value should be finite. In fact. 



(9-1) 



27r 



\q-l 2 



2(9-1)7(3-9) 



(39) 



for 1 < g < 3. The fact that, in contrast with the mean 



square displacement. 



is finite, seems to indicate 



that the constraint ( p2D is a natural one in the frame 
of anomalous-diffusion random walks Note more- 

over that Eq. (^9|) along with ( ^2| ) gives the connection 
between the Lagrange multiplier (3 and the characteristic 
length (7, 



P (X a 



(40) 



Equations (35) and ( |36| ) make clear that, as stated 
above, the maximization of entropy within Tsallis' for- 
malism does not lead to Levy distributions for the jump 
probability. Rather, a plain power-law function of the 
jump length x is obtained. Levy distributions are how- 
ever reobtained when considering the temporal evolution 
of the random walk generated by p{x). In fact, the dis- 
placement r of the walker after t time steps is given by the 
sum of the successive jumps. By virtue of the ge nerali zed 
central limit theorem (l^Jlpl discussed in Section |ll B| the 



probability distribution P{r, t) is thus given, for large t 
and 7 < 2 (i.e. q > 5/3), by a stable Levy distribution 
with Levy exponent 7. For 7 > 2, on the other hand, 
the usual form of the central limit theorem holds and the 
total displacement distribution is a Gaussian. The dy- 
namic exponent z of the random walk -which coincides 
with 7 for 7 < 2 (see Section II B)- is then 



2 ifg<5/3, 
{3-q)/{q-l) ifq>5/3. 



(41) 



This connection between the dynamic exponent and the 
Tsallis exponent q -which is illustrated by the curve d = \ 
in Figure ^- constitutes indeed the main result of the 
description of anomalous diffusion in the frame of the 
Tsallis' formulation. It shows that a close relation ex- 
ists between the properties of a Levy-flight process and 
the non-extensiveness of the involved statistics. As far as 
the underlying statistical frame differs from Boltzmann- 
Gibbs', the maximum-entropy formalism produces a ran- 
dom walk which models superdiffusion as a Levy flight. 

As in the case of the jump distribution, the total mean 
square displacement (r^) associated with Levy flights di- 
verges. On the other hand, the q-expectation value of 
is well defined for all relevant g (1 < g < 3), cf. Eq. (|3S|). 
This can be calculated taking into account the scaling 
properties of P{r,t), Eq. (pO|), and reads 
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forbidden region 
•albatkjss flight 

• DNA walk 

superdiffusion 



1 2 3 

Tsallis exponent q 

FIG. 4. The diflusion dynamic exponent z as a function of 
the TsaUis exponent q, for different spatial dimensions. The 
dots stand for some of the instances of anomalous diffusion 



discussed in Section [I A 



12 3 4 

Tsallis exponent q 

FIG. 5. Phase diagram in the [q, d)-plane, displaying the 
zones of normal and anomalous diffusion, and the forbidden 
region where the jump probability distribution cannot be nor- 
malized. The dots stand for some of the instances of anoma- 



lous diffusion discussed in Section II A 



f i:i(g)/3-H(3-9)/2 if g < 5/3, 

{r^), = <^ (42) 
y Diq)^'^^^^ ifg>5/3. 

The proportionality factor D{q) depends on q only. In- 
terpreting now the Lagrange multiplier /3 as the inverse 
of the temperature -as prescribed in the frame of Tsallis 
thermodynamics the above equation can be seen as 
a generalization of the Einstein relation (H) Equa- 
tions (Q) and imply in fact that (r^) oc /3^Tfor nor- 
mal diffusion, and Eq. (|||) is the extension of this result This connection is represented graphically in Figure ^. 



identified in a phase diagram. The connection between 
the dynamic exponent and the Tsallis exponent reads 
now 

2 if g < (4 + d)/(2 + d), 

2/(g - 1) - rf if (4 + rf)/(2 + rf) < g < (2 + d)ld. 

(45) 



to Tsallis statistics. Again, the fact that (r )g is finite for 
Levy flights suggests that Tsallis' formalism provides a 
natural frame for the statistical description of such kind 
of anomalous diffusion. 

Though the algebra is more involved than above, 
anomalous diffusion in more than one dimension can be 
straightforwardly treated in the frame of Tsallis statis- 
tics, and the main conclusions are qualitatively the same 
as for the onc-dimensional case. Maximazing the entropy 
given in ( |3l] ) in the d-dimensional space produces for- 
mally the same jump distribution as in (|33|), where the 
partition function has however to be calculated as a d- 
dimensional integral. The jump probability can be nor- 
malized if 



IV. NONLINEAR DIFFUSION AND TSALLIS 
STATISTICS 



The diffusion equation, 
dP _ 



(46) 



q < 



(43) 



and the associated mean square displacement is finite if 

i + d 



which governs the evolution of the probability P(r, t) dr 
of finding a Brownian particle in a neighborhood dr of 
point r at time t, can be generalized to take into account 
additional mechanisms acting on both the microscopic 
dynamics of the particle and the mesoscopic dynamics of 
an ensemble of such particles. A rather straightforward 
generalization introduces for instance the effect of an ex- 
ternal force field, F(r,f), acting on each particle P,^. 
This force field enters the diffusion equation as a drift 
term, namely. 



q < 



2 + d 



(44) 



Below this value, thus, the random walk generated by 
the jump probability models normal diffusion, whereas 
for (4 + d)/{2 + d) < q < {2 + d)/2 the walker performs 
superdiffusion. In Figure || these different regimes are 



dP 
'dt 



-Vr • (FP) + DVlP. 



(47) 



This equation, which combines the effect of probability 
drift -due to the force- and of probability spreading - 
due to diffusion- can be seen to govern a huge class of 
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random processes, in a generic space of states r ||3|,|31[. It 
is generally refered to as the Fokker-Planck equation. 

Further generalizations, mainly justified on a phe- 
nomenological basis, have lead to propose a nonlinear 
version of the Fokker-Planck equation P4] , namely. 



(48) 



where /i > and i' are suitable real constants. With- 
out loosing generality, one can fix /i — 1, by changing 
ptJ- p and u/ jjb — > V. In such case, Eq. (^) can be 
phenomenologically interpreted as a Fokker-Planck equa- 
tion where, ii v ^ 1^ the diffusion coefficent depends on 
the probability P{r,t). This density-dependent diffusiv- 
ity represents nonlinear effects arising, for instance, from 
interaction between the diffusing particles. Such kind 
of nonlinearities have been observed in several real pro- 
cesses, such as transport in porous media (z/ > 2) [p3[ , 
surface growth (y = 3) p^ , liquid film spreading under 
gravity — 4) and Marshak radiative heat transfer 
{1^ = 7) [|6|, among others 

The Fokker-Planck equation (^) is linear. This im- 
plies that for many forms of the force field F(r, t) the 
general exact solution can be found analytically. More- 
over, even if analytical solutions are not available or diffi- 
cult to obtain, very efficient computational methods can 
be implemented to obtain numerical solutions. In con- 
trast, exact solutions to the nonlinear equation (^) are 
particularly scarce | ]3^ , p9| , and numerical methods are 
typically subject to instabilities when dealing with non- 
linear problems. It is therefore of great interest to verify 
that, as shown in the following, Tsallis' formalism pro- 
vides a method to obtain special solutions to Eq. (^) 
for general fi and v, at least, in its one dimensional ver- 
sion and for some special forms of the drift force. This 
problem has recently been treated by Tsallis himself | |4C| ] 
and, in an alternate form, by Compte and Jou 

Consider then the one dimensional version of Eq. (^) 
for the probability density P{x,t), 



QPH 



d d^P" 



(49) 



with X = r, and with F{x) = ki — xk2. This time- 
independent, linear form of the drift force corresponds, 
in general, to a quadratic potential -i.e. an Ornstein- 
Uhlenbeck random process whereas it reduces to a 
linear potential -namely, a constant force- for k2 = 0- 
For fi = v = 1, this equation can be straightforwardly 
solved, for instance, by Fourier-Laplace transforming. A 
special solution is 



P(x,t) 



where 

Pit) _ r z(o) 



1 



exp{(3{t)[x-m?}, 



(50) 



[(1 - A)exp(-2fc2t) + A]-^ (51) 



with A = 2D/?(0)/fc2, and 

x(t) = K + [x{0) — k] cxp(— fc2i) 



(52) 



with K — ki/k2. This particular solution has the impor- 
tant property that, for t — > and /3(0) — > 00, it reduces to 
a delta-like distribution, P{x, 0) = 6[x — a;(0)]. Since Eq. 
( ^ ) is hnear for = = 1, and delta distributions can be 
used as a base for the space of initial conditions P{x, 0), a 
suitable linear combination of functions of the form ( pOj ) 
provides the solution to the linear equation for any initial 
condition in that space. In this sense, ( |50| ) gives the gen- 
eral solution to the linear Fokker-Planck equation with 
the above prescribed drift force. 

Focus now the attention on the functional form of the 
particular solution to the linear problem given in Eq. 
(^). As a function of x, P{x, t) is a Gaussian, essentially 
of the same type as (p7|). The differences are, firstly, that 
the spatial coordinate x is shifted by an amount x{t). The 
Gaussian is therefore centered around a position which 
depends on time. Secondly, the width of the Gaussian, 
which is proportional to depends also on time. 

Since the solution ( pO[ ) preserves normalization, the nor- 
malization factor 2^ is time-dependent. 

The Gaussian profile of P{x,t) in Eq. (|5^) suggests 
that this solution can be formally derived from a suitably 
extended maximum-entropy formalism, in its standard 
Boltzmann-Gibbs version. In fact, it can be shown ]4C| ] 
that such form of P{x, t) derives from the maximization 
of 



S[P] 



+ 00 



P{x,t) h\P{x,t) dx, 



with the extended constrains 



P{x,t) dx = 1, 

-00 



+00 



[x - x{t)]Pix,t) dx = 0, 



and 



+00 



[x - x{t)fP{x,t) dx 



(53) 



(54) 



(55) 



(56) 



for arbitrary x{t) and /3(t). The special forms of these 
functions that make the probability distribution satisfy 
the linear Fokker-Planck equation can be obtained by 
simply replacing P(x, t) in the equation. 

It should be by now clear that one is immediately in- 
terested at which solutions are obtained if, instead of the 
standard maximization principle, the Tsallis' formalism 
is used. Namely, take the generalized entropy 



1 



1-g 



1 



-l-CXD 



P{x,tydx 



(57) 
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and maximize it with respect to t) imposing the gen- 
erahzed constrains 



+00 



and 



[x - x{t)]P{x,tYdx = 0, 



+00 1 



(58) 



\x-x{t)YP{x,tfdx 



(59) 



This produces |40[| 



(60) 



to be compared with Eq. (|3^)- Remarkably enough, 
this form of P(a;, t) turns out to be a solution to the one- 
dimensional nonlinear Fokker-Planck equation for the lin- 
ear drift force if 



1 + j-i - ly, 



and 



m 
m 



2(1 



(61) 



(62) 



As for Levy-flight anomalous diffusion, P{x, t) can be 
normalized only if q < 3. This defines a forbidden region 
for /X > 2 + 1/ (Fig. I). 

The function Zq{t) is cxplicitely given by 



Zg{t) - Zg{0) [(1 - Ag) CM-t/^) + A< 



with 



A„ 



2iyDP{0)Zq{0Y- 



(63) 



(64) 



and r = fj,/k2{p + v)- The function x{t) is the same as 
for the linear case. Eg. (|52|). Note that the normaliza- 
tion constraint, Eq. (|54|), has not been imposed in the 
maximization of Sq[P]. In fact, the preservation of the 
norm of P(a;, t) is now not compatible with the other two 
constraints. Rather, it turns out that the integral of the 
probability density over the whole space varies with time 
according to 



+00 



P{x^ t) dx 



Ml 



+ 00 



P{x,0) dx. (65) 



This implies that the norm is conserved for all times only 
if /z = 1, or if Aq = 1 -when Zq does not depend on time. 
If Aq > 1 the norm monotonically increases for /.t > 1 
and decreases for /i < 1. If A^ < 1, on the other hand, 
the opposite behavior is observed. Moreover, for v < Q 
the norm diverges or vanishes at a finite time. Figure 
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FIG. 6. Evolution of the norm in the solutions to the non- 
linear Fokker-Planck equation for u = 1 and some values of 
fj,. These curves correspond to A, > 1. 



^ illustrates these different regimes iov v — 1 and some 
values of /i. 

The case of constant force, k2 = 0, can be analyzed as 
the suitable limit of the above solution for ^2 — >■ 0. In 
particular, taking exp(— t/r) w 1 — i/r in Eq. ( |63| ) it is 
found that 



Zq{t) - Z,(0) 



2 ^> Dmzqior-'t 



(66) 



In this limit, the width of the distribution -which is 
proportional to exhibits a well-defined power-law 

dependence on time. In fact, according to Eq. (|6^), 
l3{t)~^/^ oc t^^/i'^+i^\ This makes possible to assign a 
dynamic exponent to this kind of diffusion, given by 



z = l + -. 



(67) 



Whereas, as expected, the case = ly = 1 corresponds 
thus to normal diffusion, v/fj, > 1 corresponds to sub- 
diffusion and v/^ < 1 corresponds to superdiffusion. 
Note that ii ly/fi < —1, z < and the distribution 
width decreases with time. This seemingly unphysical 
situation |^ can however be associated with a kind of 
"negative diffusivity" , which has been observed in some 
real nonequilibrium self-organizing systems p8[ . Figure 
shows the evolution of the distribution width for /i = 1 
and various values of v. The phase diagram of Figure 
^ summarizes the various regimes obtained in different 
regions of the (v, /i)-plane. 

Equation mVj makes evident that the non-extensivity 
inherent to Tsallis statistics is related, in the frame of 
its application to the resolution of the nonlinear Fokker- 
Planck equation (^), to the nonlinearity of the equa- 
tion itself. This nonlinearity translates, at the level of 
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FIG. 7. Evolution of the width in the solutions to the non- 
linear Fokker-Planck equation for /i = 1 and some values of 
V. For the case of "negative diffusivity," /i = 1/2. 

the solutions, into anomalous properties of the involved 
transport processes. Thus, a clear connection between 
anomalous diffusion and the non-extensivity of the under- 
lying statistics arises again. It is important to point out 
that, in contrast with Eq. (po|), the solutions ( |60| ) to the 
nonlinear Fokker-Planck equation derived from Tsallis' 
formalism cannot be combined to give a general solution. 
In fact, due to the nonlinearity of Eq. (ji^), no superposi- 
tion principle holds, and ( |60|) are particular solutions for 
special initial conditions only. Nevertheless, the straight- 
forward way in which these solutions have appeared as 
an extension of the linear case along the lines of Tsal- 
lis' generalization, reinforces strongly the close relation 
between Tsallis' formalism an anomalous diffusion. 



V. CONCLUSION 

Though normal diffusion is ubiquitous in Nature, a 
large -and still growing- class of real systems is driven 
by a different kind of transport processes, namely, by 
anomalous diffusion. In view of the current importance 
of many of these systems -which range from turbulent 
flows, to disordered media and chaotic dynamics, to 
flight patterns in birds- it is of high interest having at 
hand a formulation able to place anomalous diffusion in 
a statistical-thermodynamical frame, generalizing thus 
Einstein's theory for normal diffusion. However, within 
Boltzmann-Gibbs statistics "the wonderfuU world of clus- 
ters and intermittencies and bursts that is associated 
with Levy distributions would be hidden from us if we 
depended on a maximum entropy formalism that em- 
ployed simple traditional auxiliary conditions" HJ]. It 
has been here shown that, instead, Tsallis generalized 
statistics is a strong candidate to succesfully yield such 
a formulation. 



FIG. 8. Phase diagram in the (i/, ^)-plane, displaying the 
different regimes of anomalous diffusion and norm evolution 
in the solutions to the nonlinear Fokker-Planck equation. In 
the forbidden region the probability distribution cannot be 
normalized. 

Tsallis statistics provides a natural frame for the math- 
ematical foundations of anomalous diffusion in two forms. 
In the first place, jump distributions of random-walk 
models for Levy-like superdiffusion can be straightfor- 
wardly derived from a maximum-entropy principle within 
the generalized theory. In fact, such distributions exhibit 
power-law long tails, which are an essential feature in the 
results of the theory. At once, Tsallis statistics furnishes 
an elegant explanation for the appearence of Levy distri- 
butions in other natural phenomena, as the result of the 
superposition of random variables with long-tailed dis- 
tributions. In the second place, the functional form of 
the distributions resulting from Tsallis' formalism suc- 
cessfully suggests the solution to the nonlinear Fokker- 
Planck equation, which describes both subdiffusion and 
superdiffusion. 

In order to widen the applications of Tsallis' theory to 
the statistical foundations of anomalous diffusion, further 
research should focus on some problems that still wait 
to be treated in the frame of the generalized statistics. 
For instance, it would be important to extend the deriva- 
tion of random- walk models of anomalous diffusion to the 
case of subdiffusion. As explained in Section 0, this re- 
quires introducing suitable waiting-time densities. Since 
power-law functions can fulfill this role, Tsallis statistics 
is again a natural stating point to derive such densities. 
Another extension would regard diffusion processes on 
fractals. In fact. Levy-like diffusion anomalies are the 
consequence of mechanisms driving the dynamics of the 
diffusing particles. An alternative formulation, which is 
relevant to many applications, takes into account that 
such anomalies originate rather in the complex geometry 
of the medium where particles diffuse. The connection 
between fractal geometry and Tsallis statistics has been 
identified early, and it can thus be expected that diffu- 
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sion on fractal substrates finds a satisfactory statistical- 
mechanical frame in such theory. Finally, it would be 
interesting to descend a level further in the dynamical 
bases of anomalous diffusion, and try to apply Tsallis' 
formalism to the formulation of deterministic mechani- 
cal approaches to this kind of transport. 

As a final remark it is worth mentioning that, very re- 
cently, Tsallis' formalism has been improved by redefin- 
ing the normalization of g-expectation values . This 
has solved, in a single step, two main drawbacks of the 
theory. In fact, in its original formulation, Tsallis sta- 
tistical mechanics is not invariant under energy shifts, 
and the g-expectation value of a constant depends on 
the state of the system under study. Although the cor- 
rection to the theory does not involve important changes 
in the qualitative results, it represents a major improve- 
ment from a formal viewpoint. Here, Tsallis statistics 
has been applied to anomalous diffusion in its original 
form. A relevant step forward would be to reanalyze this 
process in the frame of the corrected theory. 
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